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ABSTRACT 

The purpose of this study is to infer the effects of 
second stage engine characteristics on trajectory shape, 
weight at cutoff, and optimum control histories of the 
second stage for a two stage vehicle fulfilling circular 
orhit mission criteria. The first stage was taken to be 
an S-I booster and the second stage Is defined by an 
initial thrust to weight ratio, F/w q , and a specific Impulse, 
I S p, value. Circular orbit missions were considered at 100, 
200, and 300 km altitudes. 

By the consideration of the value of a parameter at 
the terminal point of the trajectory, we can systematically 
study the geometrical shape of the trajectory. 

Weight at cutoff is shown as a function of P/w Q for 
different I S p values, circular Injection altitudes, and 
first stage tilt programs. 

Examples of optimum control histories for the second 
stage are given. 
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SUMMARY 

The purpose of this study is to infer the effects of 
second stage engine characteristics on trajectory shape, 
weight at cutoff, and optimum control histories of the 
second stage for a two stage vehicle fulfilling circular 
orbit mission criteria. The first stage was taken to he 
an S-I booster and the second stage is defined by an 
initial thrust to weight ratio, F/w 0 , and a specific impulse, 
I S p, value. Circular orbit missions were considered at 100, 
200, and 300 km altitudes. 

By the consideration of the value of a parameter at 
the terminal point of the trajectory, we can systematically 
study the geometrical shape of the trajectory. 

Weight at cutoff is shown as a function of P/w Q for 
different Igp values, circular injection altitudes, and 
first stage tilt programs. 

Examples of optimum control histories for the second 
stage are given. 
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SECTION I. PROBLEM DESCRIPTION AND OUTLINE OP APPROACH 


We are to study the effects of second stage engine 
characteristics on trajectory shape, weight into orbit, 
and optimum control histories of a two stage vehicle ful- 
filling circular orbit mission criteria. In particular, 
the first stage is taken to be an S-I booster flying in 
two dimensions a one parameter set of angle of attack pro- 
grams in the environment of a spherical non-rotating earth 
with atmosphere. The second stage initial conditions are 
obtained by adding the earth's rotational effect at cutoff 
of the first stage under the assumption of a ninety degree 
azimuth at launch. The second stage is to fly a two dimen- 
sional calculus of variation path, minimizing the time of 
flight, in the gravitational field of a spherical earth 
free of atmosphere. In a (r, v, 5 .) coordinate system, 
the particular three first stages we are to study give 
initial conditions for the second stage of 

y 0 - 76.87 (km) y 0 = 70.37 (km) y Q = 63.21 (km) 

v 0 = 2884 (m/s) v Q = 2937 (m/s) v Q = 2991 (m/s)? 


( 1 ) 


$0 = 63.8 (deg) $o = 68.1 (deg) 


$o = 72.4 (deg) 


where y is the altitude, y = r - R OJ R 0 being the radius of 
the earth. We will refer to these first stage end conditions 
by their earth-fixed theta values, = 60 , 65 , and 70 degrees, 
respectively. 

We wish to connect these initial conditions, (l), through 
a continuous burning second stage characterized by an initial 
thrust to weight ratio, P/w 0 , and a specific impulse, I S p, 
value to circular orbit conditions at altitudes of 
y c = 100, 200, and 300 km. The variation in lift-off weight 
produced by varying the weight of the second stage in the 
form P/w 0 is ignored. The selection of a P/w 0 , Isp* and 
y c leads to a discrete minimal time calculus of variation 
second stage trajectory. The method of constructing such a 
trajectory is a numerical isolation procedure determining 
the initial values of the angle of attack and its time 
derivative, clq and d 0 . This is a numerical procedure for 
solving the two point boundary value problem. The four 
parameters of the study are taken to have the following 
restriction on their values : 
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\ 

5-j = 60 , 65 , and 70 (deg) 

.5 S F/w 0 S 2.0 I (2) 

Igp ^ 1500 (sec) 

y c = 100, 200, and 300 (km) 

J 

The following section deals with the physical shape of 
the trajectory. 


SECTION II. TRAJECTORY SHAPE 

This section deals with the geometrical shape of the 
trajectory. In order to explain the classification, let 
us consider the equations of motion for the second stage. 
In a (r, v, $•) coordinate system, 

r = v cos 

v = P/m cos a - g cos & 

& = Tp/m sin a + (g - -—) 

where P is the thrust magnitude, m is the instantaneous 
mass, and g is the magnitude of the local gravity force. 
Differentiating the first of ( 3 ) gives 


they are as follows : 

(3) 

sin $1 


= v cos 5 . - v 5- sin $ 


(4) 


Now let us examine the second stage trajectory as it 
enters the injection or cutoff point. The conditions for 
circular injection are 



$ c = 90 deg and 


g 


(5) 
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where r c is the radius of injection. Now let us evaluate 
(3) and (4) at this point common to trajectory and orbit. 
We obtain from (3) 


P 

• • - 1 - 
r c = 0 , v = — cos ot. r 
m c 


and 


P 

sin 

m v c 


(6) 


and (4) becomes 


•• 



v. 


(7) 


Combining (6) and (7) we have 


P 

m" 


sin a c . 


(8) 


Now assume we have applied the calculus of variations 
to minimize the time of flight subject to the differential 
equation constraints, (3), and with the particular values 
of y c > F/w 0 , and I S p ehoosen, the solution to the two 
point boundary value problem yields a trajectory with an 
angle of attack at cutoff, a c , of zero. Equation (8) 
states r c = 0, or that the trajectory is flat in the 
neighborhood of cutoff; here we assume that the range on 
the trajectory is a monotonic increasing function of time. 

If we assume conditions which lead to an no > 0, hence 
r c < 0, the trajectory is concave down at cutoff. In this 
case the vehicle entered the circular orbit from below. 

In a similar manner if <* c < 0 we have a concave up tra- 

jectory entering the orbit from above. 

Next, let us introduce a notation to express the general 
character of the trajectory shape: we say a trajectory is 

of class T n if it has crossed the desired free flight orbit, 
in this case a circle, exactly n times prior to cutoff. 

Since in our study all initial altitudes of the second stage 
are less than the lowest circular altitude considered (100 km), 
the T 0 , T 2 , T 4 ... trajectories inject from below and have 
“c > 0. Similarly the T a , T 3 , T 5 ... trajectories inject from 
above having <*0 < 0. Let us further denote a class dividing 
trajectory (cv c = 0) by Tij, where this trajectory divides 
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trajectories of class Ti and of class Tj (i<j). Foe example, 
class T 0 has a c > 0, and class Ti has c^c < 0, then the 
dividing trajectory would be designated by T 0 i and have 
Q'c = 0. Here, of course, we are assuming the necessary 
functional properties between points on the trajectory and 


For discussion let us fix the values of y c , $i, and 
F/w 0 and examine the family of trajectories produced by 
allowing Isp to take on different values. The study then 
proceeds by considering o-c as a function of I S p, noting the 
values of Isp for which a c = 0. These are the dividing 
trajectories. Here we should say that 

d ot c 

— ^ 0 . 

d Isp 


Figure 1 gives, as an example, this information for 

y C = 200 km, = 65 deg, and F/w 0 = .5 . It is noted 

that two such dividing trajectories are encountered in the 
range of I S p investigated, a T 01 at I S p = 160 and a T i2 
at I S p = 625. Thus, we have the following classes of 
trajectories: 


To , Isp < ' l^O 

Ti , 160 < I sp < 625 

T 2 , 625 < Isp ^ 


where the upper limit for T 2 type trajectories is not 
encountered. Figure 2 gives altitude as a function of 
range for this case. The broken lines indicate dividing 
trajectories, and the solid lines indicate typical examples 
for the classes encountered. Figures 3 through 11 give 
a c as a function of Igp for values of the parameters con- 
sistent with (2). These figures are to be compared with 
Figures 12 through 26 which illustrate the trajectory shapes. 

The case of F/w 0 = .5 yields extreme convergency problems 
in the isolation scheme for higher values of I sp : thus, some 
of the information for such cases is missing. For all values 
of F/w 0 , I sp was taken as low as practically possible. 
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In the next section we give the weight at cutoff 
information. 


SECTION III. WEIGHT AT CUTOFF 

In this section we are to investigate cutoff weight 
for values of the parameters consistent with (2). Let us 
consider the normalized cutoff weight* w c /w 0 * where wc is 
the actual cutoff weight and w 0 is the initial weight of 
the second stage* as a function of F/w 0 . This information 
is presented in Figures 27 through 47. Each figure corre- 
sponds to a different combination of the altitude of 
injection* yc, and the Igp value* while the curves of each 
figure are characterized oy a first stage* $j. 

We observe in many cases that there occurs a maximum 
of wc/w 0 with respect to F/w 0 for fixed y c * I S p* and 
values. The optimization of Wc/w Q with respect to is 
not presented due to the small range of considered. 

Although not included in this report* wc/w G as a 
function of I S p reveals a monotonic increasing behavior. 

The following section gives examples of the control 
histories . 


SECTION IV. CONTROL HISTORIES 

Here* we give typical examples of the angle of attack 
(») and chi (xT as a function of time on the trajectory. 

Chi is the angle measured from the vertical at launch to 
the vehicle longitudinal axis and is related to a as 

X = <t> + + a , 

where <p is the central angle. Figures 48* 49* and 50 give 
examples of a as a function of time and Figures 51 52* 
and 53 give the corresponding x functions. It is noted 
that for high F/w 0 and high Igp values x is nearly linear 
while this is not true for either low F/w 0 or low Igp values. 


(deg 
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y C = 200 km; *i - 65 deg; F/W 0 - .5 

ANGLE OF ATTACK AT CUTOFF (<* c ) AS FUNCTION OF SPECIFIC IMPULSE 

FIGURE 1 










Altitude (km) 
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c = 200 km; $ _[ = 65 deg; F/W 0 - . f 
GEOMETRY FOR THREE TRAJECTORY TYPES 
FIGURE 2 
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FIGURE 8 
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FIGURE 9 
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Angle at Cutoff for R c = 300 km and = 70° 


FIGURE 11 
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TRAJECTORY SHAPES FOR F/W Q = 1.0 AND ^ = 65 DEG 


FIGURE 19 
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TRAJECTORY SHAPES FOR F/W 0 = 1.5 AND = 65 DEG 

FIGURE 20 
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FIGURE 21 





TRAJECTORY SHAPES FOR F/W Q = .5 AND $j = 70 DEG 
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TRAJECTORY SHAPES FOR F/W 0 - 2.0 AND - 70 DEG 


FIGURE 26 
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FIGURE 29 
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FIGURE 46 
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R c = 200 KM = 65 DEG 

FIGURE 48 
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FIGURE 53 
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